Phase-dependent magnetoconductance fluctuations in a chaotic Josephson junction 
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Motivated by recent experiments by Den Hartog et al., we present a random-matrix theory for the 
magnetoconductance fluctuations of a chaotic quantum dot which is coupled by point contacts to 
two superconductors and one or two normal metals. There are aperiodic conductance fluctuations 
as a function of the magnetic field through the quantum dot and 27r-periodic fluctuations as a 
function of the phase difference of the superconductors. If the coupling to the superconductors 
is weak compared to the coupling to the normal metals, the (^-dependence of the conductance 
is harmonic, as observed in the experiment. In the opposite regime, the conductance becomes a 
random 27r-periodic function of 4>, in agreement with the theory of Altshuler and Spivak. The 
theoretical method employs an extension of the circular ensemble which can describe the magnetic 
field dependence of the scattering matrix. 
PACS numbers: 74.80.Fp, 05.45. +b, 74.25. Fy 
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The conductance of a mesoscopic metal shows small 
fluctuatioua of universal size /h as a function of mag- 
netic field.tl These universal conductance fluctuations are 
sample-specific, which is why a plot of conductance G 
versus magnetic field B is called a "magnetofingerprint" . 
The magnetoconductance is sample-specific because it 
depends sensitively on scattering phase shifts, and hence 
on the precise configuration of scatterers. Any agency 
which modifies phase shifts will |Hiodify the magnetocon- 
ductance. Altshuler and Spivakcl first proposed to use a 
Josephson junction for this purpose. If the metal is con- 
nected to two superconductors with a phase difference 
(f) of the order parameter, the conductance G{B,(j)) con- 
tains two types of sample-specific fluctuations: aperiodic 
fluctuations as a function of B and 27r-periodic fluctua- 
tions as a function of 0. The magnetic field should be 
sufficiently large to break time-reversal symmetry, other- 
wise the sample-specific fluctuations will be obscured by 
a much stronger B- and. (/)-dependence of the ensemble- 
averaged conductance.EfQ p 

In a recent Letter, Den Hartog et alJa reported the 
experimental observation of phase-dependent magneto- 
conductance fluctuations in a T-shaped two-dimensional 
electron gas. The horizontal arm of the T is connected 
to two superconductors, the vertical arm to a normal 
metal reservoir. The observed magnitude of the fluctua- 
tions was much smaller than e'^ /h, presumably because 
the motion in the T-junction was nearly ballistic. Larger 
fluctuations are expected if the arms of the T are closed, 
leaving only a small opening (a point contact) for elec- 
trons to enter or leave the junction. Motion in the junc- 
tion can be ballistic or diffusive, as long as it is chaotic 
the statistics of the conductance fluctuations will only 
depend on the number of modes in the point contacts 
and not on the microscopic details of the junction. 

In this paper we present a theory for phase-dependent 
magnetoconductance fluctuations in a chaotic Josephson 
junction. We distinguish two regimes, depending on the 
relative magnitude of the number of modes M and N 
in the point contacts to the superconductors and normal 



metals respectively. For M ^ N the (/>-dependence of the 
conductance is strongly anharmcuiic. This is the regime 
studied by Altshuler and Spivak.0 For M < N the oscil- 
latioig are nearly sinusoidal, as observed by Den Hartog 
et alJj The difference between the two regimes can be 
understood qualitatively in terms of interfering Feynman 
paths. In the regime M ^ N only paths with a single 
Andreev reflection contribute to the conductance. Each 
such path depends on (j) with a phase factor e^"^/^. Inter- 
ference of these paths yields a sinusoidal 0-dependence of 
the conductance. In the opposite regime M ^ N, quasi- 
particles undergo many Andreev refiections before leav- 
ing the junction. Hence higher harmonics appear, and 
the conductance becomes a random 27r-periodic function 
of 0. 
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FIG. f . Josephson junction in a three-terminal (a) and 
four-terminal (b) configuration. 

The system under consideration is shown schematically 
in Fig. It consists of a chaotic cavity in a time-reversal- 
symmetry breaking magnetic field B, which is coupled 
to two superconductors and to one or two normal met- 
als by ballistic point contacts. The superconductors (Si 
and S2) have the same voltage (defined as zero) and a 
phase difference 0. The conductance of this Josephson 
junction is measured in a three- or four-terminal config- 
uration. In the three-terminal configuration (Fig. |l|a), a 
current / flows from a normal metal Ni into the super- 
conductors. The conductance G = I/Vi is the ratio of 
/ and the voltage difference Vi between Ni and Si, S 
This corresponds to the experiment of Den Hartog et al 
In the four-terminal conflguration (Fig. Illb), a current / 
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flows from a normal metal Ni into another metal N2. The 
conductance G = I/{Vi — V2) now contains the voltage 
difference between Ni and N2. Xtis is the configuration 
studied by Altshuler and Spivak.H 

Following Ref. || we split the conductance G{B, cj)) = 
Go{B) + Grj,{B, (p) into a </>- independent background 

Go(i?)- r^GiB,cb), (1) 
Jo 

plus 27r-periodic fluctuations G^. In the absence of time- 
reversal symmetry, the ensemble average {G{B, cj))) = 
(G) is independent of B and cj). Hence (Go(i?)) = (G) 
and {G^{B, (f)) = 0. The correlator of G is 

C{5B, 6(j>) = {G{B, ^)G{B + 5B,(j) + 6(j))) - {G)\ (2) 

Fluctuations of the background conductance are de- 
scribed by the correlator of Go , 

GoidB) = (Go(B)Go(B + SB)) - {Gf 

^T'^GiSBM). (3) 

(In the second equality we used that (G^Gq) = 0.) The 
difference G^ = G — Go is the correlator of G^, 

C^SB, Scj)) = {G4B, c^)G4B + SB,c^ + 6(j))). (4) 

We compute these correlators for the three- and four- 
terminal configurations, beginning with the former. 

In the three-terminal configuration, the cavity is con- 
nected to three point contacts (Fig. |l|a). The contact to 
the normal metal has N propagating modes at the Fermi 
energy, the contacts to the superconductors have M/2 
modes each. The (N + M) x [N + M) scattering matrix 
S of the cavity is decomposed into M x M {N x N) re- 
flection matrices r (r') and N x M [M x N) transmission 
matrices t {t'), 

The conductance at zero temperature is determined by 
the igajtrix She of scattering amplitudes from electron to 
hole,B'Ll 

G = 2trs,^4, (6a) 
s^,, = -It* (1 + e'^re-"^r*y'e'^t'. (6b) 

The diagonal matrix $ has diagonal elements $„„ = 0/2 
if 1 < n < M/2 and -<^/2 if 1 -I- M/2 < n < M. We 
measure G in units of 2e^//i. 

For chaotic scattering without time-reversal symme- 
try, the matrix S is uniformly distributed in the unitary 
group.u This is the cirjCular unitary ensemble (CUE) of 
random-matrix theoryfl The CUE does not specify how 
S at different values of B is correlated. The technical 



innovation used in this work is an extension of the CUE, 
which includes the parametric dependence of the scatter- 
ing matrix on the magnetic field. The method (described 
in detail elsewhereEB) consists in replacing the magnetic 
field by a time-reversal-symmetry breaking stub (see Fig. 
H). This idea is similar in spirit to Biittiker's method pf 
modeling inelastic scattering by a phase-breaking lead.tll 
The stub contains iVstub modes. The end of the stub is 
closed, so that it conserves the number of particles with- 
out breaking phase coherence. (Biittiker's lead, in con- 
trast, is attached to a reservoir, which conserves the num- 
ber of particles by matching currents, not amplitudes, 
and therefore breaks phase coherence.) We choose our 
scattering basis such that the A'stub x ^stub reflection 
matrix rstub(B) of the stub equals the unit matrix at 
B = Q. For non-zero magnetic fields we take 

,,,,,(B)=e^^, a^^Y^^lm, (7) 

n<m 

where the matrix A is real and antisymmetric: j4„„ — 
A*^ = — Particle- number is conserved by the stub 
because rgtub is unitary, but time-reversal symmetry is 
broken, because rgtub is not symmetric if B ^ 0. In order 
to model a spatially homogeneous magnetic field, it is 
essential that -/Vstub ^ N+M. The value of A'^gtub and the 
precise choice of A are irrelevant, all results depending 
only on the single parameter a. 

The magnetic-field dependent scattering matrix S{B) 
in this model takes the form 

S{B) = Un + C/12 [1 - r3tub(S)[/22]"' r,t^biB)U2i. (8) 

The matrices Uij are the four blocks of a matrix U rep- 
resenting the scattering matrix of the cavity at i3 = 0, 
with the stub replaced by a regular lead. The distri- 
bution of U is the circular orthogonal ensemble (COE), 
which is the ensemble ofj-uniformly distributed, unitary 
and symmetric matrices.Ll The distribution of S{B) re- 
sulting from Eqs. (^ and (^) crosses over from the CQE 
for B = to the CUE for B ^ cx). One can showEl 
that it is equivalent to the distribution of scattering ixi^ 
trices following from the Pandey-Mehta HamiltoniaiiE£| 
H — Hq + iBHi [where Hq (Hi) is a real symmetric 
(antisymmetric) Gaussian distributed matrix]. 




FIG. 2. Schematic picture how the magnetic field is in- 
cluded in the scattering-matrix ensemble. A chaotic cavity 
with a spatially homogeneous magnetic field (left diagram) is 
statistically equivalent to a chaotic cavity in zero magnetic 
field (right diagram), which is coupled to a closed lead (a 
stub) having a non-symmetric reflection matrix. 
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It remains to relate the parameter a to microscopic 
properties of the cavity. We do this by computing the 
correlator Y.mn{5B) = {S,„JB)S:^^{B + SB)) from Eq. 
(^. Using the diagrammatic method of Ref. |l^ to per- 
form the average over the COE, we find (for N + M ^ 1) 



S™„ = (iV + Af)- 



[1 + (SB/B,)' 



n m, (9) 



with Be = a^^\/N + M. This correlator of scatter- 
ing matrp. elements has also been computed by other 
methods. ElrO Comparing results we can identify 



— ce^WF^^ niin(£, L)/?i(5, 



(10) 



with c a numerical coefficient of order unity depending 
on the shape of the cavity (linear dimension L, mean free 
path £, Fermi velocity vp, level spacing S). For example, 
for a disordered disk or sphere (radius L S> the coeffi- 
cient c = n/8 for the disk and 7r/15 for the sphere. 

We now proceed with the calculation of the correlator 
of the conductance. We consider broken time-reversal 
symmetry {B ^ Be) and assume that N and M are both 
3> 1. Using the method of Ref. |3| for the average over 
U, we obtain the average conductance (G) = 2NM/ {N + 
2M) and the correlator 



C{5B, 64>) = 1QKN^M^{N + M)^{N + 2M)''^ 
{N + M)2 + (iV2 + aPK) cos^{6<j)/2) 
^ [{N + M)2 - APKcos^{S(f>/2)f ' 



(11) 



where we have abbreviated K = [l + (SB/B^)'^] . Eq. 
( |lT| ) simplifies considerably in the two limiting regimes 
M <^N and M > iV. For M < iV we find 



Co{SB) = 24:{M/N)^K, 
C^{dB,6(t)) ^ 8 [M / N)"^ K COS S(j), 

whereas for M 3> -/V we have (for \S(f)\ < tt) 

-3/2 



M (5B 

'^n[-b:. 



C4SB,S<l>)^- 



^ M fSBy Af. 

IH — H 6<p 

N \BcJ 8N ^ 



(12a) 
(12b) 



(13a) 
(13b) 



The two regimes differ markedly in several respects: 

(1) The 27r-periodic conductance fluctuations are har- 
monic ii M <^ N and hi ghly a nharmonic if M ^ A^. A 
small increment d(p — ^/NjM <C 27r of the phase differ- 
ence between the superconducting contacts is sufficient 
to decorrelate the conductance if M ^ A^. 

(2) The variance of the conductance varG = Go(0) -|- 
^^(0,0) has the universal magnitude 1/2 if M A^, 
while it is reduced by a factor {8M/N)'^ if M < iV. 

(3) The variance var = C^{Q, 0) of the (/)-dependent 
conductance is larger than the variance var Go = Go (0) 
of the background conductance if M ^ (by a factor 



y^M/SN), while it is smaller if M < A^ (by a factor 
1/3). 

(4) The correlators C^{SB,0) and Co{SB) both de- 
cay as a squared Lorentzian in SB /Be if M <^ N. If 
M ^ A^, on the contrary, C^{SB, 0) decays as a squared 
Lorentzian, while Cq{SB) decays as a Lorentzian to the 
power 3/2. 

The difference between the two limiting regimes is il- 
lustrated in Fig. |[ The "sample-specific" curves in the 
upper panels were computed from Eq. (^) for a matrix 
S which was randomly drawn from the CUE. The corre- 
lators in the lower panels were computed from Eq. ([Tl|). 
The qualitative difference between M N (Fig. |^a) and 
M:^ N (Fig. |b) is clearly visible. 
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FIG. 3. Top panels: conductance minus the ensemble 
average (in units of 2e^ /h) as a function of the phase difference 
between the superconductors. Bottom panels: normalized 
correlator c{5(t>) = C(0, 5(t))/C{Q, 0), computed from Eq. (0). 
(a) is for N = 120, M = 60; (b) is for = 10, M = 160. 

We now turn to the four-terminal configuration (Fig. 
|l|b) . The two point contacts to the superconductors have 
M/2 modes each, as before; The two point contacts to the 
normal metals have A^/2 modes each. The conductance 
is given by the four-terminal generalization of Eq. (^) ,E3 



he 



(14a) 



= «/ieCjsLc», Ri! = tr SeeCjsleCi, (14b) 

)"'e**t'. (14c) 



Here (ci)m„ = 1 if m = n < A^/2 and otherwise, and 
C2 — 1 — ci. The matrix She was defined in Eq. (|6^). 
Performing the averages as before, we find (G) = N/A 
and 

C{5B, S(j)) = ^N^K[{N + M)2 + M^K cos^ (^0/2)] 

X [{N + M)^ - M^K cos^{S(j)/2)]-^ . (15) 
In the regime M <C A^ this simplifies to 
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B, 5(j)) ^ ^ {M/Nf cose 



_3_ 

32 



(16a) 
(16b) 



while in the regime M ^ N we find again Eq. ([l3| ) (with 
an extra factor of 1/16 on the r.h.s.). 

The four-terminal configuration with M 3> -/V is sinu 
ilar to the system studied by Altshuler and Spivak.El 
One basic difference is that they consider the high- 
temperature regime fcsT ^ ?i/TdwcU (with TdwcU the 
mean dwell time of a quasiparticle in the junction), while 
we assume T = (which in practice means fcsT <g; 
?i/Tdwcn)- Because of this difference in temperature 
regimes we can not make a detailed comparison with the 
results of Ref. ^. 

The features of the regime M ^ N in the three- 
terminal configuration agree qualitatively with thCpjex- 
perimental observations made by Den Hartog et al.El In 
particular, they find a nearly sinusoidal (/>-dependence 
of the conductance, with C^{B,0) being smaller than 
Co{B), while having the same i?-dependence. The mag- 
nitude of the fluctuations which they observe is much 
smaller than what we find for a point-contact coupling 
with A'l and N of comparable magnitude. This brings us 
to the prediction, that the insertion of a point contact 
in the vertical arm of the T-junction of Ref. |^ (which 
is connected to a normal metal) would have the effect of 
(1) increasing the magnitude of the magnetoconductance 
fluctuations so that it would become of order e^/h; (2) 
introducing higher harmonics in the (/)-dependence of the 
conductance. This should be a feasible experiment which 
would probe an interesting new regime. 

In conclusion, we have calculated the correlation func- 
tion of the conductance of a chaotic cavity coupled via 
point contacts to two superconductors and one or two 
normal metals, as a function of the magnetic field through 
the cavity and the phase difference between the super- 
conductors. If the superconducting point contacts domi- 
nate the conductance, the phase-dependent conductance 
fluctuations are harmonic, whereas they become highly 
anharmonic if the normal point contact limits the con- 
ductance. The harmonic regime has been observed in 
Ref. H, and we have suggested a modification of the ex- 
periment to probe the anharmonic regime as well. We 
introduced a novel technique to compute the magneto- 
conductance fluctuations, consisting in the replacement 
of the magnetic field by a time-reversal-symmetry break- 
ing stub. This extension of the circular ensemble is likely 
to be useful in other applications of random-matrix the- 
ory to mesoscopic systems. 
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